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We show that there are three possible phases for a spin- 
2 spinor Bose condensate, one more compared to the spin-1 
case. The order parameters of these phases are the sponta- 
neous magnetization and the singlet pair amplitude. Current 
estimates of scattering lengths show that all three phases have 
realizations in optically trapped alkali atoms. There is also a 
one-to-one correspondence between the structure of a spin-2 
spinor Bose condensate and that of a d-wave BCS superfluid. 

One of the recent major developments in Bose-Einstein 
Condensation (BEC) in atomic gases is the study of di- 
lute Bose gases with internal degrees of freedom. The 
first realization of such system is found in optically 
trapped ^^Na, which is a spin-1 Bose gas [Q. Recently, 
JILA has also created a "spin"-l/2 Bose gas by contin- 
ually cycling between the F — 1 and F = 2 states of 
magnetically trapped ®^Rb ||]. In the case of spin-1 Bose 
gas, the nature of the spinor condensate depends crucially 
on the magnetic interaction. In zero magnetic field, the 
spinor condensate can be either ferromagnetic or "polar" , 
which has very different properties |^ |Q] . 

Generally, only atoms in the low lying hyperfine multi- 
plet are confined in the optical trap. Those in the higher 
hyperfine multiplet will leave the trap by spin-flip scatter- 
ing. In case of ^^Na and *^Rb, their hyperfine multiplets 
{F — 2 and F = 1) are regular, i.e. the higher spin state 
{F = 2) has higher energy. Since spin-flip scattering is 
strong in ^'^Na, it may be difficult to produce a spin-2 
Bose gas in this system. On the other hand, ^^Rb has 
much weaker spin-flip scattering and is a candidate for 
optically trapped spin-2 Bose gas. In the case of ^^Rb, 
the lower multiplet has spin F — 2. It also has a negative 
s-wave scattering length in zero field. Should the current 
effort to Bose condense ®^Rb in magnetic traps be suc- 
cessful, it is conceivable that an = 2 spinor condensate 
can be trapped optically in low fields, if the three parti- 
cle losses when the field is reduced through the Feshbach 
resonance is not too large. 

In this paper, we study the ground state structure of a 
spin-2 Bose gas within the single condensate approxima- 
tion. In the case of spin-1 Bose gas, it has been realized 
recently that the ground state can be "fragmented" , (i.e 
containing more than one condensate) ^ Q. Despite 
this fact, the phase diagram for single spinor conden- 
sates remains highly valuable and in fact gives the best 
agreement with experiments so far This is because 
the spin-1 fragmented state is delicate with respect to 
spin-nonconserving perturbations, which will drive the 
system toward a single condensate state. For these rea- 
sons, we shall first focus on the ground states of single 
spinor condensates. We shall consider only linear Zee- 



man effect, which is already much more subtle than the 
spin-1 case. The actual fragmented structures as well as 
quadratic Zeeman effects will be discussed elsewhere. 

Because of the increase in spin value, spin-2 Bose gas 
has one more interaction parameter than that of spin-1 
Bose gas. As a result, there are three possible phases in 
zero magnetic field (instead of two in the spin-1 case). 
These phases are characterized by a pair of order param- 
eters (f) and O describing the ferromagnetic order and 
the formation of singlet pairs respectively. The order pa- 
rameter O is absent in the spin-1 case. The first two 
phases, which are characterised by (|(f)| = 2, 8 = 0) and 
((f) — 0, |0| = 1) in zero field, are referred to as ferro- 
magnetic and polar phases, respectively. They are the 
analogs of the corresponding phases in the spin-1 case. 
The third phase, ((f) = 0,0 = 0) is a non-magnetic 
but degenerate set of states which has no spin-1 analog. 
They will be referred to as the "cyclic" states because of 
their close analog to the d-wave BCS superfluids which 
we shall discuss at the end. In finite fields (along z), 
both (fz) and O are non-zero for all three phases. From 
the current estimates of scattering lengths, we find that 
the spin-2 Bosons of ^^Rb and ^^Na will have a "polar" 
ground state, whereas those of *^Rb and ^^Rb will be 
cyclic and ferromagnetic respectively. All these phases 
can be distinguished by their very different density pro- 
files. § 

Low-energy Hamiltonian. The effective low-energy 
Hamiltonian of a spin-/ Bose gas was derived earlier |^ 
with particle interactions of the form V(ri — r2)=5(ri — 

Y1'f=o9f'^f^ 9f ~ 47r?i^a_F/M, where M is the mass 
of the atom, Vf is the projection operator which projects 
the pair 1 and 2 into a total hyperfine spin F state, 
and ap the s-wave scattering length in the total spin 
F channel. Symmetry implies that only even F terms 
appear in V. For spin / = 2 Bosons, we have V = 

tVi + g2V2+9oVo. Using the fact that Vo+V2+Vi = 1 
and fi • f2 = (F^ - if - i^)/2, V2 and Va can be ex- 
pressed in terms of fi •f2 and Vq. The resulting expression 
for the interaction is: 

V(ri - V2) = S{vi - V2) {a + I3h ■ h + 57^o) , (1) 



where a = \{-ig2 H 
K50-54)- 1(52-54) 
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is then 



354), fi = -7(52 - 54) and 7 = 
The second-quantized Hamiltonian 
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+57V'+V'a' (2a; 2a'|00) (00|2&; 2byj,^,, 



(2) 



(2/3(/,)-p)/,C-AC + 27eiC = 0, 



(6) 



where (00|2&;26') is the Clebsh-Gordon coefficient for 
combining two spin-2 particles with = b and b' 
into a spin singlet |0, 0), '0a(r) (a — 2, . . . — 2) an- 
nihilates a Bose at point r with spin a, and /C = 

/ '^'^ ( W^V'l^V'a + VtrapV'lV'a) is the sum of kinetic en- 
ergy and trap energy Vtrap- The Po term represents the 
linear Zeeman shift. It also includes a Lagrange multi- 
plier pi to constrain the total spin so that in an external 
field B along z its expression is po = pi + gusBh If 
the ground state \G) is un- fragmented, the field operator 
-00 (r) in eq.(||) becomes a c- number 'l'a(r) = (V'a(r)) — 
a/ n(r)Ca(r), where n is the density and the C is a nor- 
malized spinor, C2Ca — 1- The ground state energy then 
becomes 

{H)g ^K- j npoih) + \Jn^a + + j\Qf) 

(3) 

where K = (/C)g, (f> = C*^fC, the superscript 
stands for transpose, and e = v^(00|2&;26')C6Cb'- More 
exphcitly, 9 = Ca^afcCfc = 2C2C-2 - 2CiC-i + Co> where 
Aab = ^a-i-6,o(— 1)°- Note that 6 represents a singlet 
pair of identical spin-2 particles and is therefore invari- 
ant under any rotation U = e^^^" (where c is the ro- 
tational angle), i.e. C^^C = C'^U'^AUC- This implies 
("^ {f'^ A + Ai)C, = for arbitrary C, which in turn implies 
that Af is antisymmetric, and hence d^'^Ai^ — 0. 

It is useful to note that Q is the scalar product of a 
state C with its time-reversed state ( = A(* , i.e. 9 = C^C- 
That C is the time-reversed state can be seen from the 
fact that C^{C=C^A^fAC=-C^i^C = -C^fC- When C 
is equal to its time-reversed partner up to a phase factor, 
(i.e. A(* = a^, \a\ = 1), we say that C has no broken 
time-reversal symmetry. In this case, |8| = 1 Time- 
reversal symmetry is broken if A(^* ^ a( for any \a\ = 1. 

Determination of Spin Structure. The first step to un- 
derstand the spin structure is to consider the homoge- 
neous case, where the spin configuration ( is controlled 
by the energy function 



5(C) = /3(f>^ + 7|ep-p(/.), (p 



2po, 



(4) 



As we shall see, once the spin structure of the homoge- 
nous system is determined, the resulting phase diagram 
will provide a quick though qualitative determination of 
the actual structure. It is easy to see from cq.(^ that 
the ground state magnetization must be aligned with the 
external field (i.e. along z), implying (/+) = 0, where 
/+ = fx +ify Eq. 



then becomes 



f-/3(/.)' + 7|0|'-p(/.)- 



(5) 



Minimizing eq.(|^), we have the following Euler-Lagrange 
equations 



where A is the Lagrange multiplier for the normalization 

By contracting eq.(^ with (and with C"^^: respec- 
tively) to the right and using the antisymmetric property 
of fA, we obtain: 



2/3(/,)2 + 27|ep-p(/,) -A = 

e(A - 27) = 0. 



(7) 

(8) 



Eq.(||) leads to the following cases: (a) 8^0, and (b) 
8 = 0. Case (a), which implies A = 2-7, will be referred 
to as the poZar phase. In case (b), eq.(||)-(0) yield 



(2/3(/,)-p)(/,-(/,))C = 0, 



(9) 



which further divides into : ih.l){fz — {fz))C = {fer- 
romagnetic phases) and (b.2) p = 2l3{fz) {cyclic). It is 
obvious that the ground state is degenerate under gauge 
transformation and spin rotation along z. The family of 
degeneracy is larger when p = Q because of the full ro- 
tational symmetry. Without loss of generality, we only 
consider p > 0. 

Polar phases : Setting A = 27 in eq.(|^) and denoting 
1 = AC. we find 



(2/3(/,)/, - pf, - 27)C + 278C - 
(-2/3(/,)/, +pA - 27)C + 278*C = 
[(2/3(/,) - pffl - 472 + 472|8p]C - 



(10) 

(11) 
(12) 



Eq. (|1 If) is obtained by multiplying ( [lOD with A to the 
right and taking the complex conjugate. Eq. (|2|) shows 
C, is an eigenstate of with possible eigenvalues = 
4, 1,0, denoted as P, PI, and PO respectively : 



PI : 



PO : C 



1 



P 



o«"2 .111 P. p*"-2 4 - - 

2 + 8/3-27'^'"'"''' V2 8/3-27 



(0,0,1,0,0), 



(13) 



where are arbitrary phases. At p = 0, 8 = 1. 

Ferromagnetic phases. Case (b.l) 8 = and (/^ — 
{fz))C = leads to two ferromagnetic phases, 

F: = (1,0,0,0,0), F': = (0, 1, 0, 0, 0) (14) 

with energies £ = 4/3 — 2p and E = [3 — p respectively. 

Cyclic phases. Case (b.2) corresponds to 8 = and 

p = 2/3(/z). From eq.(||), we see that the states satisfying 

these two conditions are degenerate with energy £c — 
2 

— -^g. Explicitly, this degenerate family is specified by 
the conditions 
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2C2*Ci + Veer Co + V6Co*C-i + 2CiC 



(/+) 

e = 2C2C-2-2CiC-i + Co -0 

p = 2/3(2|C2p-2|C-2p + |Ci|'- 

IC2P + IC1P + IC0P + IC-1P 

One example of such spinor is 



IC- 



2 = (15) 
(16) 
(17) 
(18) 



(19) 



where </> is an arbitrary phase. 

Phase diagram in zero field. Table || summarizes the 
above results. When p = 0, all polar states are degen- 
erate. By searching for the lowest energy state in Table 
I for given scattering lengths (hence given /3 and 7) , we 
obtain the phase diagram in (7, /3)-plane as shown in the 
inset of Fig.l. Among the two ferromagnetic phases only 
the state F is realized. F' is never the lowest energy state 
for all (7, /?). The three phases P, F and C are separated 
by the boundaries /3 = (F-C), 7 (C-P) and 4/3 = 7 
(P-F). 

The phases P, C, and F are characterized by the "or- 
der parameter" (|e| = 1, (/,) = 0), (6 = 0, (/,) = 0), 
and (6 = 0, (/z) — 2) respectively. Moreover, time- 
reversal symmetry is broken for the cyclic and the fer- 
romagnetic states but not for the polar states, where C 
and AC* are related by a phase factor, (see eq.jl^)) [||. 
Since the pair (6, (fz)) undergoes discontinuous changes 
from one phase to another, the transition between differ- 
ent phases at p = are all first order. 

It is also useful to display the phase diagram in terms of 
the differences in scattering lengths — and 02 — 04 
as shown in Fig.l. The regions occupied by the three 
phases are : 

P : ao - 04 < 0, ||a2 - a^] < ^\ao - 04]. 

F : (a2 - 04) > 0, i(ao - 04) + |(a2 - a^) > 0. 

C : (02 — 04) < 0, ^(ao — 04) — |(a2 — 04) > 0. 

Based on the current estimates of scattering lengths (in 

a.u.) by J. Burke and C. Greene 



Spin-2 species 


ao 


a2 


04 


^^Na 
s^Rb 
S5Rb 

83Rb 


34.9 ± 1.0 
89.4 ±3.0 

-445.0l^°g 
83.0 ±3 


45.8 ± 1.1 
94.5 ±3.0 

82.0 ±3 


64.5 ± 1.3 
106.0 ±4.0 

-420.0t}^° 
81.0 ±3 



we note from Fig.l that all the three phases can be real- 
ized in the above alkali (and in fact in Rb) isotopes. The 
error bars in the estimates of scattering lengths, however, 
introduce uncertainties in the predictions of ground state 
structures. The case of ^^Rb is particularly unclear, for 
it barely resides in the polar region while the error bar 
covers both polar and the the cyclic states. The fact 
that all these realizations are close to the phase bound- 
aries means that other physical effects such as gradient 



energy and magnetic field gradients will be important 
in determining the spatial structures of condensate, for 
they will compete with the energy differences between 
different phases. 

Phase diagram in finite field. When p ^ 0, the phase 
diagram depends on the sign of the Heisenberg interac- 
tion p. The phase diagrams for /3 > and /3 < are 
shown in Fig. 2(a) and 2(b) respectively. In all cases, only 
the polar state P is realized as PI and PO are never the 
lowest energy states. When /3 > 0, (Fig. (2a)), the ex- 
ternal field p polarizes the non-magnetic polar and cyclic 
states. Since 9 assumes different values in P and C, the 
boundary between them is first order. As p increases, 
both P and C gradually cross over to the ferromagnetic 
state. The second order phase boundary at which this 
occurs is p = 4/3 (p = 4/3 — 7) for the C (P) state. 
When /3 < 0, (Fig. (2b)), the only non-magnetic state is 
P, which will also be polarized by the external field p as 
in the /3 > case. The phase boundary between the polar 
state and the ferromagnetic state is again a second order 
line, given hy p = — (4|/3| + 7). Of course, time-reversal 
symmetry is broken for all phases when p — 0. 

Fig. 2(a) and (b) also provide a quick (qualitative) de- 
termination of the spatial spin structure of non-magnetic 
phases in optical traps. Since p is inversely proportional 
to the density n, (see eq.(^)), it generates a vertical tra- 
jectory in the phase diagram as one moves from the cen- 
ter of the trap to the surface of the atom cloud. Such 
trajectories are shown in Fig(2a) and (2b). The start- 
ing points A, B, and C indicate different kinds of spin 
states at the trap center. The fact that p diverges at the 
surface of the cloud means all non-magnetic states will 
develop an outer ferromagnetic layer whenever p ^ Q. 
The above construction does not apply to the ferromag- 
netic condensate as its spin structure is not affected by 
p. Numerical calculation including the kinetic energy is 
needed to determine the spatial spin configuration with 
fixed total magnetization. 

Because of the different "order parameters" 9 and 
(/z), these three phases have very different spinor struc- 
tures. By measuring the density in each spin component, 
they can be easily distinguished from each other. 

Finally, we note that there is a one-to-one correspond- 
ing between the structure of a spin-2 spinor Bose con- 
densate and that of a d-wave BCS superfiuids. The lat- 
ter is known to be characterized by an order param- 
eter which is a traceless 3x3 symmetry matrix By , 
{B = B^ , TiB = 0). To see this correspondence, we 
note from the property of spherical harmonics that the 
sum P(k) = /c^^m5^2m(k) is a homogenous polynomial of 
degree 2 in k (-P(k) = Bijkikj), and satisfies the Laplace 
equation V|P(k) = 0. These conditions guarantee that 
_B is a traceless symmetric matrix. We can therefore as- 
sociate with each spinor a traceless symmetric matrix 

^ / ^fci%C)n^2m(k)- Rewriting the energy in tcrms 
of B, one finds that eq.(Q) in zero field {p — 0) reduces to 
the free energy of a d-wave superfiuid. The exact mini- 



3 



mization of this problem was solved by N.D. Mermin . 
Our zero field results are in agreement with his exact so- 
lution. Because of this correspondence, we have named 
the spinor phases in zero field according to the features of 
the d-wave solutions It is also clear from the above 
discussion that the structure of a spin-S* spinor conden- 
sate {S — integer) corresponds to that of a singlet BCS 
superfluid with orbital angular momentum S. 

We would like to thank Eric Cornell and Carl Wie- 
man for discussions on their experiments on ^^Rb, and 
to Jim Burke and Chris Greene for estimates of scatter- 
ing lengths. This work was supported by NASA grant 
NAG8-1441 and NSF grants DMR-9808125 and DMR- 
9708274. 
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TABLE I. Possible phases in finite field 
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Figure Captions Fig.l. Phase diagram in zero field. 
The ferromagnetic, polar, cyclic phases are denoted as 
F, P, and C resp. The phase boundaries P — F and 
P — C are straight lines with slopes —7/10 and 7/10 resp. 
The boundary F — C is 02 — 04 = 0. All boundaries are 
first order lines. The locations of various alkali isotopes 
on this diagram are the symbols : o =^'^Na, x =^^Rb, 
A =®^Rb, • =®'^Rb. Inset: Zero-field phase diagram in 
the (7,/?) plane. 

Fig. 2. Spinor phase diagrams in finite field for (a) 
P > and (b) (3 < 0. The thick and dashed lines are 
first and second order phase boundaries resp. For an 
atom cloud with a non-magnetic state (A, B, or C) at 
the center, as one moves from the center to the surface, a 
vertical line is generated in the phase diagram, showing 
that the outer layer is always ferromagnetic. 
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